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Abstract: We report a bound on the relative frequency fluctuations in 
nonlinear second harmonic generation. A 1064 nm Nd:YAG laser is used 
to read out the phase of a Mach-Zehnder interferometer while PPKTP, a 
nonlinear crystal, is placed in each arm to generate second harmonic light. 
By comparing the arm length difference of the Mach Zehnder as read out by 
the fundamental 1064 nm light, and its second harmonic at 532 nm, we can 
bound the excess frequency noise introduced in the harmonic generation 
process. We report an amplitude spectral density of frequency noise with 
total RMS frequency deviation of 3mHz and a minimum value of 20 
/iHz/Hz 1 / 2 over 250 seconds with a measurement bandwidth of 128 Hz, 
corresponding to an Allan deviation of 10 -19 at 20 seconds. 
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1. Introduction 

Harmonic generation of optical fields in nonlinear optical crystals can be modeled with clas- 
sical coupled wave equations 1 1 ]. Define the fundamental and second harmonic fields as E\ = 
S\ exp(27liV\t) and E2 = ^exp^TT/v^), respectively. In ideal second harmonic generation, 
V"2 = 2v\. While this frequency ratio has been measured very precisely ((V2/V1) = 7 x 10 -19 
[2]), we present a lower bound to the frequency noise, or the spectrum of the time dependent 
quantity V2 — 2V\. 

Several experiments at the forefront of precision metrology and frequency standards use 
harmonic generation in their experiments. These include iodine stabilized Nd:YAG lasers (3, 
optical frequency combs (UE1, measurement of optical frequency ratios |2||6), and precision 
atomic spectroscopy |7|. Many of these experiments provide bounds to any excess frequency 
noise which might be found in the second harmonic generation (SHG) process (e.g. from ther- 
modynamic fluctuations in the crystal temperature). As experiments like these push towards 
lower noise levels, the fundamental noise sources in second harmonic generation may become 
a relevant noise source. 

Another field which utilizes precision SHG is interferometric detection of gravitational 
waves. These interferometers require measurements of mirror displacements at the level of 
10 -20 m/ V^Hz 1 8]. To achieve such high precision, these interferometers use seismically 
isolated suspended mirrors which make up multiple coupled optical cavities with narrow 
line widths. It has proven challenging to control these cavities and bring them simultaneously 
into resonance. In order to simplify this problem, auxiliary fields will be generated by SHG and 
used to form independent signals for two of the cavities. This technique is being developed us- 
ing some small scale prototypes OHO). To understand the limits of such a scheme, we need to 



bound the frequency noise introduced in the second harmonic generation process. In our case, 
we use a 1064 nm carrier (Nd:YAG), and generate the 532 nm doubled light with PPKTP crys- 
tals. In order to meet the stability requirements for this technique, the excess frequency noise 
must be less than 70 mHz RMS from 10 mHz to 30 Hz 0. In the future, optical configurations 
using more than one wavelength inside the interferometer to beat standard quantum noise lim- 
its of the detection [11] may be used. Fundamental noise limits in harmonic generation may be 
relevant there. 

2. Experimental Setup 



LEGEND 
M1 - HR 1064 with PZT 
M2- HR 1064 
M3.M4-HR 532/1064 
M5,M6-low loss HR 1064 
M7.M8 - HR 532 
BS1 -50:50 1064 
BS2- 50:50 532/1064 
G1,G2 - PD Balancing gains 
DAQ - Data Acquisition System 
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Fig. 1. Experimental layout - A dual wavelength Mach-Zehnder Interferometer. Some read- 
out mirrors for harmonic separation are not shown. See text for a more complete descrip- 
tion. 

Figure [T] shows the experimental setup used to measure the uncorrelated frequency noise 
between the fundamental and the second harmonic. The beam from a 2 W non-planar ring laser 
at 1064 nm is passed through a spatial filter cavity (l2\ . The spatially filtered beam then enters 
the Mach-Zehnder through a hole in an acoustic enclosure. PPKTP crystals (Raicol Crystals) 
inside temperature stabilized ovens are placed in each arm, with lenses added to mode match to 
the each crystal. Two additional dichroic mirrors (M3, M4) are placed in one arm for alignment. 
The beams recombine on a dichroic 50:50 beamsplitter and low loss HR mirrors (T < 10 ppm at 
1064 nm) are used to separate the fundamental from the second harmonic. Commercial dichroic 
mirrors (HR532, AR1064) were used to further separate the 532 nm from the 1064 nm light. For 
the 532 nm and 1064 nm detection, Si (Hamamatsu 1223) and InGaAs (GPD 2000) photodiodes 
were used, respectively. The Mach-Zehnder arm lengths were adjusted to be mid fringe for 
both 532 and 1064 nm simultaneously. The interferometer is locked to this point by applying 
feedback from the 1064 nm PD difference signal to the PZT on Ml (see Fig.[T]). We drove the 
PZT with a small signal at 100 Hz to monitor the calibration (rad/V) at both wavelengths, as any 
differential frequency noise between the two wavelengths would cause the 532 nm calibration 
to drift. The calibration of the Mach Zehnder was confirmed by sweeping through multiple 
fringes. The 1064 nm and 532 nm Mach-Zehnder length signals were read out with two pairs 
of balanced homodyne detectors. The sum and difference signals were digitized at 256 Hz and 



processed further offline. 

In the non-depleted pump approximation, with imperfect phase matching, the phase relation 
of the two fields in each arm before the recombining beamsplitter is: 

0532 = 201064 - n/2 - ML/2 , (1) 

where 01064 and #532 describe the relative phases of the fundamental and the second harmonic, 
respectively, L is the length of the doubling crystal, and Ak is the phase mismatch, the parameter 
normally used to describe the efficiency of SHG. In theory, Ak can be arbitrarily small (limited 
in practice only by the ability to stabilize the temperature of the nonlinear crystal). 

Using superscripts to differentiate between the two arms, the difference in the second har- 
monic phase is thus 

0532 " 0532 = 2 ( U>64 " ^) ~ " A* 8 ) L/2 . (2) 

3. Results 
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Fig. 2. Example time series of the Mach-Zehnder output. The RMS phase noise over the 
one hour period was 6 mrad RMS at 532 nm and 3 mrad RMS at 1064 nm 



The phase difference (50io64(O = #1064(0 — #fo64(0) was suppressed by the servo, which 
had a unity gain frequency of ~ 10 Hz. Typical values of 50io64(O anc * 50532(0 are shown in 
Fig-i 

The frequency noise amplitude spectral density of the Mach Zehnder is shown in Fig. [3] 
with known noise sources. <5vio64 and 5V532 are shown in red and green, respectively, where 
8 v = 86/271. Above 10 Hz the phase noise is dominated by a forest of mechanical resonances 
on the optical table which show up strongly in both <5vio64 and 5V532. We used the average 
transfer function over the measurement time, H(/) = (5v532(/)/5vio64(/))» t0 subtract the 
noise which is coherent between the two 5 vs. This subtracted level is shown as the black trace 
in Fig. [3] The black trace bounds any noise source which causes frequency noise between the 
fundamental and the harmonic and is uncorrected between the two ovens, such as thermody- 
namic fluctuations in the crystals, or any temperature fluctuation not common to both ovens. 
Some common mode effects such as intensity dependent phase shifts and temperature fluctua- 
tions (and thus phase matching fluctuations) are suppressed by the experimental setup, so these 
technical noise sources will not be visible in these measurements. However, a pessimistic esti- 
mate of the temperature noise coupling assuming no common mode rejection shown in Fig. [3] 




10 2 10 1 10° io 1 io 2 10 2 10 1 10° io 1 io 2 

Frequency (Hz) Frequency (Hz) 



Fig. 3. On the left we see the full noise budget of the experiment: 0f O64 — 0^)64 ( m rec ^ 
(which is the in-loop signal of the Mach-Zehnder), 6^ 2 ~ Q532 ( m § reen ) (which contains 
excess phase noise), and the subtraction residual (in black). The total RMS frequency fluc- 
tuation in the measurement band is 3 mHz. On the right we have a comparison of this work 
with previous bounds. These bounds show up as frequency noise [13], timing error |5], and 
phase noise (4). 

is below the measured excess frequency noise. In addition, the intensity of the fundamental was 
only stabilized to a level of 2 x 10 -6 at 1 Hz, where it can be lowered to the 10 -8 at 1 Hz with 
current techniques, so unless the rejection of this effect was more than 46 dB, we can safely 
ignore it. It is highly improbable that the intensity to frequency noise coupling would be above 
the level shown in Fig. [3] The excess noise below 10 Hz was found to be correlated with air 
currents on the table, and would be reduced by moving the setup into a vacuum chamber. The 
total RMS excess frequency noise of the black trace is 3 mHz RMS in the 10 mHz to 128 Hz 
band. 

4. Noise Sources 

In addition to the usual technical noise sources, it is worth considering whether there is a more 
fundamental limit to the relative phase between the fundamental and the harmonic. A rough 
estimate of thermal noise from thermoelastic (Zener) damping was obtained by directly ap- 
plying fill the Fluctuation-Dissipation Theorem. We treat the crystal as an 0.8 mm radius 
cylinder, and follow the calculation done by Heinert et al. 1 15 ]. This yields a spectral density of 
Ak taking into account both thermorefractive and thermoelastic fluctuations. When expressed as 
frequency fluctuations, it is well approximated by 2.5/(1 + 500/ -7 / 8 ) /iHz/ \/\\z above 10 Hz. 
Below that, it must be flat or continue to decrease, or else the RMS temperature integral would 
diverge. Practically speaking, in our band of interest, the temperature fluctuations of the ovens 
is at least 2 orders of magnitude greater than these fundamental thermodynamic temperature 
fluctuation. In the future, when researchers seek to make frequency comparisons at better than 
the 10 -21 level, these thermal noises will have to be calculated with more accuracy. 



5. Previous Bounds on Excess Frequency Noise 



We have examined the results from a number of precision experiments which involve SHG. 
Although these experiments are likely limited by other noise sources, these data can be used to 
place upper bounds on the possible noise generated by SHG. In Fig. [3] we compare frequency 
noise fT3lL timing errors | 5 ], and phase noise El; in Fig.|4]we compare Allan deviations 0[3l[5] 
[6] [T6l [T7j HS1 - While the comparison in Fig.[3]is straightforward, some caution should be taken 
interpreting Fig. [4] Our Allan deviation at the 0.1 s time scale is heavily influenced by the high 
frequency noise in the measurement (10-128 Hz). Since we low pass the signal to acquire data 
at 256 Hz, we reject noise which would make the Allan deviation increase at all time scales. It 
should also be emphasized that we only measure relative frequency fluctuations, and that there 
are some common mode noise sources which the experiment is insensitive to. See lTT9l for more 
information on Allan deviations and phase noise. 




Fig. 4. The Allan deviations from this experiment (green and black) were obtained from the 
power spectrum as described in |20|. Bounds from previous work are shown for reference. 



6. Conclusions 

In conclusion, we have demonstrated that the RMS frequency fluctuations added from uncor- 
rected mechanisms between two SHG crystals is less than 3mHz at time scales over 10 ms. 
The obvious correlated mechanisms (temperature and intensity noise coupling) are likely in- 
significant compared to this as discussed in section [3] This is low enough to not limit the lock 
acquisition O [TOl scheme for Advanced LIGO and other gravitational- wave detectors. Addi- 
tionally, we have shown that there is no excess noise process at a level which is of interest 
to those doing precision atomic spectroscopy |7], using frequency combs to transfer optical 
harmonics (HO EH and other tests of fundamental physics. 
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